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Abstract. In this article we develop geometric versions of the classical Langevin equa- 
tion on regular submanifolds in euclidean space in an easy, natural way and combine them 
with a bunch of applications. The equations are formulated as Stratonovich stochastic 
differential equations on manifolds. The first version of the geometric Langevin equation 
has already been detected before, see |LRS12| for a different derivation. We propose an 
additional extension of the models, the geometric Langevin equations with velocity of 
constant absolute value. The latters are seemingly new and provide a galaxy of new, 
beautiful and powerful mathematical models. Up to the authors best knowledge there 
are not many mathematical papers available dealing with geometric Langevin processes. 
We connect the first version of the geometric Langevin equation via proving that its 
generator coincides with the generalized Langevin operator proposed in |Sol95) . |j0r78| 
or [KolOO] . All our studies are strongly motivated by industrial applications in modeling 
the fiber lay-down dynamics in the production process of nonwovens. We light up the 
geometry occuring in these models and show up the connection with the spherical veloc- 
ity version of the geometric Langevin process. Moreover, as a main point, we construct 
new smooth industrial relevant three-dimensional fiber lay-down models involving the 
spherical Langevin process. Finally, relations to a class of self-propelled interacting par- 
ticle systems with roosting force are presented and further applications of the geometric 
Langevin equations are given. 



1. Introduction 

This article is about geometric Langevin equations on regular submanifolds of euclidean 
space and its applications. As starting point, let us consider the classical Langevin equation 

dCt = (^tdt (LI) 
dwt = -Xujt dt - V^>(6)dt + adWt 

in M^'', d G N. For simphcity <I> G C°°(M'^) and A, cr G (0,cxd). Besides its crucial applica- 
tions in statistical physics we regard this model as some kind of universal kinetic equation 
for the evolution of a random particle described by position and velocity coordinates mov- 
ing in under the influence of an external force V<I>. 

Motivated by our research and development of so called fiber lay-down models we are 



in need to derive specific manifold- valued analogues of Equation (1.1) as well as spherical 
velocity models out of the latter. Herein, fiber lay-down processes arise in the production 
process of nonwovens and the expression is used for the description of the forms generated 
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by the stochastic lay-down of flexible fibers onto a moving conveyor belt. The under- 
standing and mathematical simulation of such fiber webs is of great industrial interest, 
see |KMW09] and references therein. A simplified stochastic model simulating a virtual 
fiber in a fast and efficient way is developed in |GKMW07] and describes the lay-down of 
a single fiber as a curve in M?. More realistic three-dimensional models are derived first 
in |KMW12"a] . Consider |GKMS12] for their mathematical analysis. 

As explained in fG KMSl"2] the two- and three-dimensional model can at once be for- 
mulated with the help of the manifold-valued Stratonovich stochastic differential equation 
(SDE) with state space x E<^-\ d G N, d > 2, given by 

d6 = ^tdt (1.2) 
dut = -{I -ut® oJt) V^>(^t) di + a[I -ujt® ojt) o dWt. 

Here S*^^^ is the unit sphere in and x0y = xy^ . All basic notations as well as a short 
introduction to the concept of Stratonovich SDEs on manifolds is presented in Section |2} 
In Section [3] we start lighting up the geometry occuring behind this basic fiber lay-down 
model. We give a new interpretation and will see that this equation is just the natural 



analogue of the classical Langevin Equation (1.1 ) having spherical velocities. 



The stochastics occuring in the fiber Equation (1.2) are modeled with the help of a 
Brownian motion on S"'^^. Compared to realistic lay down-processes this gives to rough 
paths and one aims to have a model with smoother trajectories. Clearly, in one may 
simply replace a standard Brownian motion through a classical Langevin equation with 
potential $ = to obtain a smoother stochastic process. In the spherical case, it is not a 
priori clear how to get a smoother version of some Brownian motion and how to formulate 
a spherical Langevin equation. 

So this motivates to search for natural manifold- valued versions of the classical Langevin 
equation. This is done in Section |4] with the help of a simple transformation strategy 
that gives the general manifold- valued Stratonovich SDE, called the geometric Langevin 
equation^ in a surprisingly easy and natural way. The transformation strategy reads as 
follows. Via observing that the geometric Langevin equation must have state space equal 
to the phase space of some manifold M, i.e., the tangent bundle TM, we first formulate the 
Langevin equation in its most natural way on the tangent bundle of M/27rZ. This space is 
diffeomorphic to TS^. As consequence, we get a natural equation for the Langevin process 
on TS^ which can immediately translated to TS'^ and afterwards finally even easily to TM. 
Here M may be any regular manifold in of dimension d as introduced in Section [2] 
The equation with state space TM then reads 

dit = ujtdi (1.3) 
dujt = -Xujt dt - F{^t,(^t) dt - gradM ^>(^t) dt + aUM[^t] o dWt- 

Here gradfj^ denotes the (tangential) gradient on M, Hm[C] the orthogonal projection onto 
the tangent space at the point ^ G M and F((^,oj) is a suitable forcing term keeping the 
trajectories on TM. For details and notations, see Section [2] and Section [4j This equation 
has also been detected before by Lelievre, Rousset and Stoltz, see |LRS12| Eq. (3.3)] or 
|LRS10| Sec. 3.3]. Nevertheless, our strategy to obtain the natural generalization of the 
classical Langevin equation completely differs from the strategy given in [LRS12] and 
moreover, exactly the same arguments afterwards serve as starting point for developing 
the above mentioned smoother version of the basic fiber lay-down model. We also strongly 
believe that the underlying strategy is helpful for every applied mathematician who aims to 
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derive similiar manifold-valued stochastic kinetic equations. Therefore, we should present 
the approach in full detail. 

Up to the authors best knowledge there are apparently not many mathematical papers 
available dealing with the construction and analysis of Langevin type equations in its 
general geometric form, see Section [5} We found the surprisingly and seemingly unknown 
paper of Soloveitchik, see |Sol95] . in which the author constructs the Langevin process 
(called Ornstein-Uhlenbeck process therein) with external potential on smooth, compact 
and connected Riemannian manifolds. First, Soloveitchik introduces a natural generaliza- 
tion of the classical Langevin generator and constructs afterwards its diffusion process. 
In this section we connect our geometric Langevin equation with Soloveitchiks approach 



via proving that the generator L associated to (1.3) coincides with the one proposed by 



Soloveitchik, i.e., for some local coordinate system [q^ , . . . ,q'^, v^, . . . , v'^) on TM we obtain 
L=y { r-' v v Xv^ ] + -a > 

Having in mind our applications, such local generator representations then directly lead to 



convenient numerical simulation methods for equation (1.3) via formulating so called local 



SDEs, see Remark 5.2, Moreover, let us already remark the article [J0r78i by J0rgensen 
in which the same operator appears without forcing term as well as the article |KolOOj in 
which Kolokoltsov even discusses more general versions of the operator from above. 

In view of our applications we restrict attention afterwards to the spherical situation and 
discuss the spherical Langevin procss on TS^ seperately. This includes explicit simulation 
formulas and is done in Section |6l 

In Section [7] we apply the previous machinery to develop the smooth fiber lay-down 
model. The driving stochastics therein is now given by the desired spherical Langevin 
process. The general equation has state space x TS'^^^, d G N, d > 2, and reads 

dit = 0Jtdi (1.4) 
dwt = -(/ -iot® wt)V$(^t) dt + /xt dt 

d^t = • V^{it) dt - A/itdt - l^tpwt dt + a{I -ujt®^t)° ^Wt- 

The physical relevant cases are d = 2 and d = 3. This model is developed in coopera- 
tion with our colleagues Klar, Maringer and Wegener, i.e., the authors from |KMW12b] . 
The equation was first detected in its full form including all necessary terms by the sec- 
ond named author of the underlying article. This important fact is not mentioned in 
|KMW12b] . Again we provide simple local coordinate expressions guaranteeing convenient 



numerical simulations of Equation (1.4), see Proposition 7.2 A detailed comparision be- 
tween the basic and the smooth fiber lay-down model and its industrial applications will 
be discussed in a forthcoming research paper of the authors from |KMW12b] together with 
the authors of the present article. 

In Section [T] we present further applications to some spherical velocity models arising in 
self-propelled interacting systems describing the collective behaviour of swarms of animals, 
see |CKMT10] . They have a similiar mathematical structure as the two-dimensional fiber 
lay-down equation. Analogously as in the fiber lay-down scenario they can now directly 
be translated into the three-dimensional case and moreover, while substituting a spherical 
Brownian motion through a spherical Langevin process, new smooth versions of the models 
can be formulated. 
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Motivated by the basic fiber lay-down model, we derive in Section [8] a natural version 
of the geometric Langevin equation moving with velocity of constant absolute value. This 
equation is then called the geometric Langevin equation with spherical velocity. It lives on 
S^M, the spherical tangent bundle of M of radius r > 0, and is given by 

d6 = a;tdt (1.5) 
dwi = -F{iuu:t) dt - (n§,,[a;t]nM[6])v$(et) dt + a {li^A^t] Hm^]) ° dWt. 

Up to the authors best knowledge, this equation seems to be new and has nowhere proposed 
before. This equation contains in case M = M*^ the basic fiber lay-down model as special 
case. Moreover, there are now a whole galaxy of new beautiful and powerful mathematical 
models available. 

The reader can see the beauty of the spherical velocity version of the geometric Langevin 
equation in Section [9] where we discuss specific examples. One example deals again with 
the spherical situation. Therein, the spherical velocity version of the geometric Langevin 
process on S^-M, M = S^, is constructed and a special coordinate representation is derived. 
Altogether, this gives a further smooth spherical process which can again be used to 
construct a new smooth three-dimensional fiber lay-down model. We expect that the 
latter serves as real alternative to all existing three-dimensional lay-down models. An 
additional research article on this application is intended. 

We remark that the discussion about the geometry of the basic fiber lay-down in Sec- 
tion |3] serves partly as motivation. The reader who is not interested in this application 
may directly switch to Section [4] Furthermore, in Section [2] we present the necessary 
mathematical background and introduce notations that are used until the end of this 
article. 

Finally, we believe that this article is interesting from both a theoretical and applied 
point of view. Sometimes, pure and applied mathematics are seemingly disjoint. With 
the help of this article, we also hope to build a small bridge connecting some real world 
applications with the more theoretical area of SDEs on manifolds. 

2. The setup and notations 

Before starting with the geometric derivation and discussion of Langevin type equations 
on general submanifolds in the euclidean space, we introduce some notations and recall 
well-known facts. We stay detailed in order to be self-contained. First some general 
notations: C°°(M"), n E N, denotes the set of all infinitely often differentiable functions 
/ : M" — )• M. The index c means compact support. V (or V^.) always denotes the 
usual gradient operator (with respect to the variable x) as column vector and the 
Hessian matrix in the euclidean space. | • | is the standard euclidean norm. The standard 
euclidean scalar product is simply denoted by • or also by (•, Oeuc- Superscript T denotes 
the transpose and rank the usual rank of some matrix. The expression smooth means that 
the underlying object is of class C°° . I is the identity matrix. Partial derivatives with 
respect to some variable x are denoted as usual by ^ or for short by dx- Convention: 
Any vector x E is always understood as column vector. And the notation (x, y) for 
x E M"', y E M"^ is understood as 

{x,y) = 

Until the end of this article we follow the notations and language introduced in the 
underlying section without further mention this again. 



j E M"+^ 
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2.1. Submanifolds in the euclidean space. In the following, let d,k £N and define 
N := d + k. We fix some c/-dimensional (regular) submanifold M C M-^ of the form 



{eeM^ I =0, =0}, rankJ^(/i,...,/fe) =A:, 



(2.6) 



Here G C°^(M^^) and := J5(/i, 

matrix of /i, . . . , /fc at the point ^ G M, i.e., 



,fk) 



denotes the Jacobian 



Jdh^ ■■■■:fk) 



N ^ ' 



dfk 



CG 



In case fc = 1, we write H := M. In this situation we introduce the following (orientated) 
unit normal vector field n on H as 

Coming back to the general situation, the tangent space T^M at the point ^ G M, 
embedded in R'^ , is now given by 

r^M = G I J^w = O} = G I a; • Vfj{0 = 0, j = 1, . . . , A;} 

We identify each element uj from the embedded tangent space T^M C M^, £, G M, with 
its associated derivation lo ■ V(^) from the algebraic tangent space, in notation 

io = oj- V(0- 

Then Lo{f) = u ■ V/(0, f eC°°{U),U C M open with ^ G C/, is understood as co ■ Vg{^) 
where g G C°°{R^) is a function which extends / locally in some open neighbourhood of 
M at ^. Recall that oj ■ V5(^) is independent of each such extension g for /. 

Now the tangent bundle of M, i.e., TM = JJ^^iyj T^M, is given as embedded manifold in 
M^^ simply by 

TM = {{^,uj) G M^^ I fjiO =0,u- VfjiO =0, j = l,...,k}. 

So TM is again a (regular) submanifold of M?^ of dimension 2d since with hj{^, u) := fj{C) 
for j = 1, . . . ,k, and with hj{^, lo) := u ■ V fj{Cj for j = + 1, . . . , 2A;, we have for the 
Jacobian 



J{hl,...Mk)- 



Here V'^fjUJ is understood as the mapping M^-^ 9 (C)*^) '-^ V^/j(^)a; G M-'^ for all 
J = 1 . . . , A;. Consequently, we get 



rankJ(^^;^)(^i,...,^2fe) = 2A;, (^,w) G 
Furthermore, we define the spherical tangent bundle bundle S^M of radius r as 



= r^} , r > 0. 



In case r = 1 we shortly write SM and call it the unit tangent bundle. The function 
^2fe+i(C)^) '■= Y^f=i ~ ''^ serves as an additional defining function for SrM. Then it is 

easy to sec that §,.M is a (regular) submanifold of of dimension 2(i — 1. TM and S^M 
will later on serve as the right state spaces for generalized Langevin type equations. 
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Now choose a smooth mapping TM 3 {^,u}) i— t- ^(^,0;) G M . Then ^ is a smooth 
vector field on TM, i.e., it holds additionally ^(^,0;) G T(^,a;)™ for any (^,a;) G TM, if 
and only if 

(Y^)-^ = 0, (^^1") .^ = 0,^ = 1,.. .,A:. (2.7) 

Moreover, A is even a smooth vector field on S.f.M if additionally to (2.7) we have 



(^^^ ■ A{C,oj) = 0, (e,a;) G S,.M. (2.5 



Further recall that for some local coordinate chart {U,q) of M with q : U ^ q(U), 
q{U) C M'^, the vectors 

dr dr 

form a basis of T^M, ^ G M. Here r is defined as r := on g(C/) and s G q{U) is chosen 
such that r(s) = ^. Let q^ , j = 1, . . . ,d, denote the coordinate functions of q. Moreover, 
the vector field from q{U) corresponds under the local diffeomorphism q to the vector 

field ^ on U, j = 1, . . . ,d. More precisely, let ^ be defined by 

1-/(0 :=4/({)-|-/W.i = l.....rf. 

where / G C°°{q{U)) and / G C°°([/) are related by / o r = /. Here ^ is the usual 
algebraic derivation induced by the coordinate functions. Hence = ^ since 



In the same way as in (2.9) one computes corresponding vector fields under some diffeo- 
morphism between two manifolds. 

The corresponding local coordinate chart {TU, Tq) of TM is now determined through 
TU = {(C, uj) £TM \ i£U} and 

(Tq) : TU ^ q{U) x (C, u) ^{q\^),..., q^{0,v\C, u^),---, w)) 

where the are uniquely determined by the relation uj = J2'j=i ^"'(C)'^) §f:{0- other 
words, we have the diffeomorphism 

/ d 

ds 



{Tq)-^ ■.q{U)xM.'' ^TU, (s, k) ^ (e(s), c^(s, k)) := | r(s), ^ -^(s) | (2.10) 



where s = (si,...,Srf) G ^(f/) and k = (ki,...,k^) G M'^. Furthermore, M inherits 
canonically the structure of a Riemannian manifold obtained from the underlying euclidean 
space M^. Thus the Riemannian metric is determined in the coordinate chart {U,q) via 
9ii ~ §s^ ' §s~ for i, j = 1, . . . , d. The inverse matrix associated to (gij) is denoted as usual 
by (g*-'). The Christoffel symbols are given by the formula 



2 ^ y dsm dsn dSj 
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Through using the previous relation for the gij one immediately verifies 

dr d'^T 
dsi dsnds 



Km = T.9''7^^-7^n^^ ^,n,m = l,...,d. (2.11) 



Next we introduce the orthogonal projection from onto T^M, ^ G M, as 



Um[C]=I- [J^Ji) J5. (2.12) 

By the rank condition (2.6), note that ^J^ really exists and is smooth as con- 

sequence of the powerful implicit function theorem. In case M = H, (2.12) reduces to 
He = / — n n. Here x := xy^, x,y ^ R^. Finally, the (tangential) gradient of some 
/ G C°°(M) reads 

gradM /(O = nwie] V/(0, e e M. (2.13) 

The latter is again understood as applied to some smooth local extension g (in M^) of / 
around ^. This is indeed well-defined. 

Remark 2.1. In order to include itself into this framwork, we regard M"^ as canonically 
embedded in M'^+^ via 

M"' = G R'^+^ I id+i = 0}. (2.14) 

Then note that S^M, M := M"*, can canonically be identified with x S^"^ with S^'^ the 
sphere of radius r in M'^. 

Finally, we need some basics about Stratonovich stochastic differential equations (ab- 
breviated by SDEs) on manifolds discussed next. 

2.2. Stratonovich SDEs on manifolds. We shortly discuss the concept of manifold- 
valued Stratonovich SDEs. Excellent mathematical treatments on this subject can e.g. be 
found in |Hsu02] . |HT94] or |IW89j . So let X be a general (abstract) smooth manifold 
and let {Q-,F,f,{Ft}t>o) be a standard filtered probability space equipped with an r- 
dimensional standard {J^f}t>o-Brownian motion W = {Wt}t>Q- Let Vq, Vi,...,Vr be 
smooth vector fields on X and let x G X. A solution X = {Xt}t>o of the Stratonovich 
stochastic differential equation 

r 

AXt = Vo {Xt) dt + ^ Vj {Xt) o dW^'^ (2.15) 
i=i 

with initial condition Xq = x is any {J^j}t>o- adapted, continuous process on X (the 
one point-compactification of X) having A as a trap such that the following is satisfied: 
Xq = X holds P-a.s. and for every / G C^(X) the process {f{Xt)}t>o satisfies the M-valued 
Stratonovich SDE 



df{Xt) = {Vof){Xt) dt + Y.{V,f){Xs) o dW}' 

i=i 

Any such X has generator L : C°°(X) ^ C°°(X), L = Vo + 5 Ei=i V|, and is a L-diffusion 
process, see [IIsu02, Sec. 1.3]. Such L-diffusion processes are weakly unique in the sense 
that they induce a unique L-diffusion measure on the path space, see |Hsu021 Sec. 1.3]. 
Moreover, the Stratonovich solution concept behaves in a natural way under state space 
transformations: If X is diffeomorphic to some Y with diffeomorphsim 99 : X — )■ Y, then 
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Y = ip{X) solves the associated Stratonovich SDE on Y in which the vector fields Vj from 
X are replaced through their push-forward vector fields under ip, see |Hsu021 Prop. 1.2.4]. 
We call the SDEs on X and Y equivalent. 

Finally, if X C W" is a regular submanifold as defined in the previous section with embed- 



ded vector fields Vj then SDE (2.15 ) can equivalently be viewed as some usual Stratonovich 
SDE in M" in which the vector fields Vj are arbitrarily extended to smooth vector fields on 



(denoted again by Vj). This means that the solution to (2.15 ) is obtained by solving the 
extended equation in M". The solution to the extended equation then stays on X provided 
that the initial value lies on X, see |Hsu02| Prop. 1.2.8]. In this case we introduce some 
notations. Define the matrix valued mapping p ^ V{p) = (Vi(p), . . . , Vr{p)), p G IR", and 
set 



VodWt ■.= '£{Vjf){Xs)odW^ 



t ■ 



3. Starting point: The geometry of fiber lay-down and a new modeling 

point of view 

As already mentioned in the introduction, this section serves partly as motivation. The 
reader who is probably not familiar with the fiber lay-down model may skip this section 
in first reading. Nevertheless, due to didactic reasons, we start in this section with the 
discussion of the geometry occuring behind the so called (basic) fiber lay-down equations. 
In particular, we give a new geometric modeling of the latters. We recall that all of these 
equations are simplified stochastic models simulating a virtual fiber web as observed and 
needed for optimization in the production process of nonwovens. For further motivation 
we refer to the description from the introduction and to |KMW09| . 

Several mathematical questions arise out of the fiber lay-down geometry dealing essen- 
tially with questions about geometric Langevin type equations and serve as motivation for 
all our further studies in the rest of the underlying article. As described in the introduc- 
tion, the aim of developing a smoother version of the basic fiber lay-down model yields 
us to the study of manifold- valued versions of the classical Langevin equation. Moreover, 
we will see that the basic fiber lay-down model, developed originally in |GKMW07| and 
generalized in |KMW12a] and |GKMS12] to higher dimensions, is just a special realization 
of a geometric equation which will later on be called the geometric Langevin equation with 
spherical velocity, see Section [8} 

But first, we start by revisiting the original modeling of fiber lay-down. 

3.1. The arc-length modeling of fiber lay-down revisited. The two-dimensional 
fiber lay-down equation in case of a non-moving conveyor belt is originally introduced in 
|GKMW07j . The lay-down of a single fiber is modeled therein as a curve ^ in M^. The basic 
assumption and starting point of the authors is that the curve is assumed to be arc-length 
parametrized. Consequently, one has \dt(,t\ = 1 where t has just the interpretation of 
arc-length. With the ansatz dt^t = (cos(af), sin(at))^, the authors of |GKMW07] propose 
a stochastic equation in a in order to get a realistic fiber lay-down behaviour given by 
the arc-length parametrized curve ^j, t > 0. Furthermore, in all upcoming articles and 
models, see e.g. |BGK-f08| , |KMWQ9) . |KMW12a] and |KMW12bj . the interpretation of t 
being the arc-length is still valid. Consequently, considering this point of view, the models 
do not have a priori a physical interpretation. 

From a physical point of view, the evolution equation for the fiber lay-down process 
must be obtained in a natural way with t really being time. Henceforth, let us propose 



GEOMETRIC LANGEVIN EQUATIONS AND APPLICATIONS 



9 



an alternative point of view which is even fundamental apart from the fiber lay-down 
discussion, see Section [8| 

3.2. The basic fiber lay-down equation: A second view by view and a new 
modeling approach. Therefore, we consider the following scenario in the production 
process of nonwovens (see |BGK-|-08| ). We also assume the conveyor belt to be stationary. 
The spinning speed of a single fiber, i.e., the amount of (fiber) material coming out of a 
single nozzle, is assumed to be a constant equal to Vg G (0, cxd). Here the fiber may be seen 
as an infinitesimal thin one-dimensional, elastic, slender object. Chosing the time unit 
appropriately we may assume without loss of generality that Vs = 1- The fiber lay-down 
process now draws a curve G M'^, t > the time index, on the stationary belt. Here 
either d = 2 or d = 3 are treating the physical relevant scenarios. Due to an turbulent flow 
in the deposition region near the conveyor belt, the behaviour of the curve is of stochastic 
nature. By the previous it is now reasonable to assume that the lay-down speed of the 
fiber (material) again coincides with its spinning speed Vs = 1- This includes that we 
assume the fiber to be inextensible in the non-moving conveyor belt case. So, if a;^ denotes 
the attached velocity vector at the lay-down curve at time t > 0, it must hold |a;t| = 1. 
Thus the tuple {^t,^t) lives on the unit tangent bundle of M*^ for all t > 0. The latter is 
equal to M*^ x where S*^"^ denotes the unit sphere in M"^. Moreover, in reality one 
has an ergodic behaviour of the fiber lay-down process (in case of a stationary belt) and 
observes that the distribution of converges towards a unique stationary state as t — t- oo. 



Now consider a particle which is governed by the classical Langevin equation (1.1). Then 



its density in phase space M^*^ fulfills the associated Fokker-Planck evolution equation. A 



stationary solution to the latter is e *(^)e . Thus one expects that the stationary 



distribution of the particle is given up to normalization by e *(^)e dxdv. Because of 



this we regard the Langevin SDE (1.1 ) as some kind of universal stochastic kinetic model 
for an ergodic process described by position and velocity coordinates. 

By the previous discussion we propose the most simpliest ansatz for the ergodic stochas- 
tic fiber lay-down process from above in form of an usual Langevin SDE (1.1) in M^'^ for 
(^t,a;f), t > 0. The potential <I> G C°°(M'^) then has to be determined later on. The further 
restriction \ujt\ = 1, i > 0, suggests to do a suitable spherical projection of the velocity 
coordinates u) in this equation. Intuitively, one has to project an infinitesimal step of the 



wt-component in (1.1) onto the sphere. This is analogous to the construction of spherical 
Brownian motion, see |RW871 Ch. 5, page 183]. This projection of course has to be be 
done with the Stratonovich operation which obeys the natural transformation rules from 
ordinary calculus, see e.g. |Hsu02] . In other words, the describing Stratonovich equation 
for the easiest case of the fiber lay-down process with state space M'^ x S'^^^ simply reads 

d(,t = ujt dt 

dwt = -(/ - wt (8) ujt) {Xujt + V5^>(^t)) dt + a{I -ut^ ut) o dWt- 
Here W is standard d-dimensional Brownian motion and A, a are both nonnegative con- 



stants. Furthermore, recall the notation (/ — w (8> (^) o dWt introduced in Section 2.2 The 
solution to this equation stays on X 

1 provided that the initial value lies on the 
latter manifold, see |Hsu02[ Prop. 1.2.8]. Also note that (/ - a; u;)uj = for a; G S'^~^ 
Consequently, the equation reduces to 

d^t = ujtdt (3.16) 

dujt = -(/ -iot® u)t)V^<^{(,t) dt + a {I -uJt<^ ut) o dWt 
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with state space M"' x S'^^^. This is exactly the basic two-dimensional (for d = 2) and 
three-dimensional (for d = 3) fiber lay-down equation, see |GKMSl"2] for details and 
explicit simulation formulas. For the original derivation consider |KMW12a] . But now 
t is really the time and the equation comes out in a completely easy and natural way, 
without any calculation. As shown in [GKM ST2] . a stationary distribution to the latter 
equation is given by e^^'^^^^^d^, S, d € N, d > 2, where S denotes the Riemannian 
volume measure on S"^^-*^ and d^ the usual Lebesgue measure in M.'^. Here $ is usually 
chosen as ^>(^) = ui^f + o"2^2 '^'^s^i^ '^i > 0, ensuring that the lay-down curve comes back 
to its reference point determined by the nozzle. 

Finally, we remark that the curve 7 which describes the lay-down in case of a non- 
stationary conveyor belt with speed Vb G [0, 1] moving in direction e, |e| = 1, may now also 
naturally be modeled as = + Vbt e, t > 0. This is a slightly different and alternative 
modeling suggestion as considered by the arc-length point of view in |BGK-|-08] . 

Summarizing, apart from the fiber lay-down scenario, the discussion altogether shows 
that Equation (3.16) may be interpreted as some kind of a universal stochastic kinetic 
model for an ergodic process described by position and velocity coordinates moving in M"^ 
with velocity of constant absolute value. This seems to be fundamental and interesting by 
itself and in Section [8] we use and generalize the previous projection strategy once more. 



4. Derivation of the geometric Langevin process on regular submanifolds 

AND ITS StRATONOVICH SDE 

As starting point, let us consider the Langevin type equation 

dxt = vtdt (4.17) 
dvt = -A vt dt - Vx'i'ixt) dt + ao dWt 

with state space M^'^, (i G N. Here ^' G C°°(M°') is a suitable potential function, A and 
cr are nonnegative constants and W denotes a d-dimensional standard Brownian motion. 
Clearly, adWt = aodWt- As mentioned in the introduction we imaginize this equation as 
some kind of general kinetic model and still call it the Langevin equation. This process is 
placed at time t > at the position xt G M*^ and has an attached velocity vector vt G M'^, 
thus its state space is interpreted as the tangent bundle of W^. In the following we derive 
a natural manifold-valued analogue of the Langevin equation. For the same reasons as 
before this process must live on the tangent bundle TM. Here the manifold M C can 
be chosen in its general form as defined in Section [2j The final equation is then called the 
geometric Langevin equation. 

We explicitly remark again that our geometric Langevin equation has already been 
detected in its form before, see |LRS12| Eq. (3.3)]. Nevertheless, the strategies for deriving 
the equation are different. The original motivation for us to derive a geometric Langevin 
equation arised during the development of a three-dimensional smooth fiber lay-down 
model, see the motivation from the introduction of this article. The strategy to obtain 
the desired geometric Langevin equation is summarized and explained in the introduction. 
We make use of the latter again in Section [7j 

Furthermore, we remark that the transformation strategy below is similiar to the one in 
|GKMS12j . In the latter article we have given a mathematical precise geometric derivation 
of the (original) three-dimensional fiber lay-down model from |KMW12a] . Consequently, 
original ideas for such a strategy can be found in |KMW12a] and have been the starting 
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point for all our differentialgeometric considerations. Now let us start first with the for- 
mulation of the Langevin equation on M/27rZ x M. We recall that all necessary statements 
concerning Stratonovich SDEs on manifolds are summarized in Section 2.2 see |Hsu02| 
for details. 



4.1. The Langevin equation on R/27rZ x M. Now we formulate the Langevin equation 
(4.17) in its most natural way on the tangent bundle of M/27rZ, i.e., on M/27rZ x M. Any 
point p G M/27rZ x M is written in the form p = (a, v). We define the abstract Stratonovich 
Langevin equation on M/27rZ x M as 



dXt = AoiXt) dt + Ai{Xt) o dWt 



(4.18) 



where W = {Wt}t>o is a standard one-dimensional Brownian motion and the vector fields 
-^1 are defined by 

d 



V — (p) 



d d 



Here $ G C°°( 
tion 



i/27rZ). Let us justify the latter equation. Consider the canonical projec- 



P : 



V27rZ, x^a = [x]. 



Now take any solution Yt = (xt, Vf), t >0, to the Langevin equation (4.17) in in which 
the potential ^' is chosen as ^' := $ o P. Let / G C°°(M/27rZ x M) and let g £ C°°{R'^) be 
given by g{x,v) := f{P{x),v). Define X hy Xt = {at,vt), at = P{xt), t>0. We have 



dg df 
9{Yt) = f{Xt), T^{x,v) = — (Q,t;), {x,v) G 



P(x). 



dx ' da 

Thus together with the Stratonovich transformation rule (see e.g. |Hsu02j ) we obtain 



<if{Xt) = dg{Yt) 



vt ^{Yt) dt - {Xvt + d^^ixt)) ^{Yt) dt + a^{Yt) o dWt 
{Aof){Xt)dt + {Aif){Xt)odWt. 



In other words, X solves equation (4.18). This justifies Equation (4.18) on the tangent 
bundle of R/27rZ in a natural way. 

4.2. The Langevin equation on TS^ and TS'^. Now ]R/27rZ x M is diffeomorph to TS^, 
where S"*^ is the unit sphere embedded in M^. Consequently, the Langevin equation on TS^ 
can be derived via calculating the equivalent SDE to Equation (4.18) on TS^. But first 
note 



{(e,a;)GM^| I^P = 1, u; ■ ^ = o} 



Hence V is a smooth vector field on TS^ if and only if 



V(e,a;) 



0, 







(4.19) 



holds for all (^,w) G TS^. The previously mentioned diffeomorphism from R/27rZ x R to 
T§^ is now given as 

ip{a) \ ( \ _ fcos{a) 



- sin(a) 
cos (a) 
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which is induced by the diffeomorphism Lp : M/27rZ — ^ S-*^. Analogously as in Section |2j 
see Equation (2.9), we get 



da 



-V ip 



d_ 
dv 



Here the tilde denotes the push- forward vector fields on TS^. Define ^ G C°°(S^) via 



$ := $ 



-1 



. Then 

d<^> 

da 



-6 

6 



Thus ^P'^ = gradgi<&(^) since T^S^ = spanj^-*"}. Furthermore, v'^ip = Thus 



^0 = 



Iwp^- Aa;-gradsi$(0, 



^1 ^ 



Ao and Ai are defined in the previous section. So the equivalent Stratonovich SDE to 
(4.18) for the Langevin equation on TS^ reads 

dXt = MXt) dt + MXt) o dWt (4.20) 

where is a standard one-dimensional Brownian motion. Observe that -^i really 
satisfy the conditions in (4.19). In order to obtain a natural translation of the Langevin 
equation to TS'*, S*^ the unit sphere in M'^+-^, we slightly change the stochastic part in 
(4.20). Therefore, we replace o dWt through 

2 

{i-c<^o°dWt = Y^ 



(I 







where ej, j = 1,2, denotes the j-th unit vector in and W now denotes a standard 
two-dimensional Brownian motion. The resulting Stratonovich SDE on TS^ then differs 
from (4.20), but it is easy to see that its generator is exactly the same as the one from 
(4.20). In particular, any solution to the new modified equation on TS^ then coincides 
weakly with any solution to (4.20). Now we may write the resulting equation, directly 

(4.21) 



translated to TS'^, simply as some Stratonovich SDE in K2(d+i) in 

the form 

d^t = oJt dt 

dujt = -Xutdt - Iwtl^^dt - gradgd ^>(6) dt f7n§d[^t] o dWt- 



The solution to the latter equation stays on TS"' provided that the initial value lies on 
the manifold. Here is a standard (d -|- l)-dimensional Brownian motion and recall 
nsd[C] =1 — ^^Cj^^ S*^- By abuse of notation, the potential is written again without 
tilde. 

4.3. The geometric Langevin equation on TM. Let M C be as in Section [2] Hav- 
ing in mind the Langevin equation on TS"' from (4.21 ) we propose a natural ansatz for the 
geometric Langevin equation with state space TM as the following abstract Stratonovich 
SDE 



N 



dXt = Ni{Xt) dt + N2{Xt) dt + N^{Xt) di + Y, MjiXt) o dWi 



U) 



(4.22) 
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where is a standard A^-dimensional Brownian motion and the vector fields Ni and ^Aj 
are defined on T 



as 



w 



Mo 





-\oj 





-gradM $(0 



a 



n. 







(4.23) 



Here e,-, j 



1, 



term F : T. 
TM. By condition (2.7) this means 

72 



,N, is the j-th unit vector in M.^ , $ G C°°(M), and the smooth forcing 
^ has to be chosen such that the vector vield A/i is really tangential to 



So the tangential condition (2.7) holds by setting 

VVi(0^^ 



a; 



F{c,u) :=jnj^jjy 



(4.24) 



Note that A/2, As and all Mj indeed satisfy the tangential condition (2.7). As described 



in Section 2.2 the solution to the abstract geometric Langevin Equation (4.22 ) is obtained 
by solving the following Stratonovich SDE in M?^ 

d£,t = ujt dt 

dujt = -A uJt dt - F{^t, uJt) dt- grad^ ^{^t) dt + a Umi^t] 

whose solution stays on 
occuring vector fields from 



(4.25) 

[6] o dWt 

I provided that the initial value lies on TM. Of course, the 
in (4.25 ) have to be understood as arbitrary smooth exten- 



sions of the vector fields from TM in (4.22 ). Thus our modeling of the geometric Langevin 



equation is complete. At this point we mention again [LRS12j or [LRSIO . Sec. 3.3] for 
an alternative, different derivation of the geometric Langevin process in the context of 
constrained Langevin dynamics. 



Remark 4.1. Furthermore, note that in case 



or 



we get back (4.21) or 



(4.17) respectively. In the second case recall the identification of M made in (2.14). Of 



course, we neglect afterwards the redundant {d + l)-component of ^ and u in the resulting 
Langevin equation in the M'^-case. Finally, note that the embedded SDE (4.22) writes as 
Ito-SDE in in exactly the same form. 

5. The Langevin generator 

Let us cite the mathematical articles dealing with generalized geometric Langevin pro- 
cesses. Besides the already mentioned article |LRS12] (or |LRS10| Sec. 3.3] respectively), 
consider the book of Gliklikh, see |Gli97] , for the construction and discussion of Langevin 
type equations with external forcing field arising in geometric mechanics. In the latter, the 
Ito SDE approach (involving Ito bundles) on general Riemannian manifolds is used. Fur- 
thermore, consider |j0r78| where the Langevin process on the tangent bundle of smooth 
manifolds is constructed via the Gangolli-McKean injection scheme. Therein, J0rgensen 
provides a general expression in local coordinates for the generator associated with the 
Langevin process. A forcing field (or external potential) is not included there. And as 
described in the introduction, Soloveitchik, see |Sol95] . also considers a construction ap- 
proach via defining a suitable generalized Langevin (or called Ornstein-Uhlenbeck therein) 
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generator L. This generator extends the one from J0rgensen and includes the additional 
potential term. Finally, consider the book of Kolokoltsov, see |KolOO] . in which the au- 
thor constructs a curvilinear Ornstein-Uhlenbeck process on the cotangent bundle of some 
compact d-dimensional manifold via introducing a suitable operator extending the one 
from Soloveitchik, see |KolOO| Ch. 4, Sec. 1]. Some of the above mentioned articles are 
seemingly not known to each other, see the references therein. 

In the following, we calculate the generator associated to our previously defined geomet- 
ric Langevin process and show that it coincides with the one introduced by Soloveitchik. 
Such a local representation is not contained in |LRS12j , or |LRS10j respectively. This con- 
nects both approaches and moreover, yields a desirable way to formulate (4.22) in local 



coordinate form, see Remark 5.2 



2.1 



Proposition 5.1. LetWi he a d- dimensional regular submanifold o/M^ as in Section 
assume that ^ E C°°(M) and let (JJ,q) be a local coordinate chart o/M. The Kolmogorov 
operator L : C~(TM) ^ C°°{T 

TM is given in the corresponding coordinate chart (TU,Tq), Tq 
on TM as 



associated to the geometric Langevin equation (|4.22) on 

.4\ 



,q'^,v^ 



dqi 



d 



9- 



dq^ dv^ 



d 



1 



dv'^dv^ 



Proof. Step 1: First remember the local diffeomorphism {Tq) '■ q{U) x 



TU map- 



ping {s,k) to {T{s),J2j f^jdjT^s)), see (2.10). Here djT 



Under {Tq) the smooth 



vector fields 

asi 



d 



on q{U) X M'^ correspond to the smooth vector fields 



d_ 
dsi 



d_ 
dq^ 



d 



d_ 



from TU . This follows in the same way as in (2.9). 



Step 2: The Langevin generator L is given on C°°(TM) as 







(5.26) 



see Equation (4.23). We write each (^,a;) G TU again in the form (t(s), KjdjT{s)) with 
{s, k) G q{U) X M*^ as introduced in Section 



2.1 



Let us calculate the vector fields Mj and 
Aii in the local coordinate system {TU, Tq) . We start with A/i. By the previous, there 
exists uniquely determined aj(^, w), 6j(^, w) G M such that 



-F{(,u;: 



For notational convenience we omit the argument {(,,u)) in the following. By (5.26) and 



the uniqueness of the representation for uj we conclude that aj = kj must hold for all 



j = 1, . . . ,d. Now by taking the scalar product with respect to some {^q.^J both sides 
of the latter equation and using again ( 5.26[ ) we get 

Kn Km diT ■ dnmT + E ^" ^"'^ ' = 
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Here we have used the fact F ■ diT = which fohows since J^v = holds for ah v £ T^M 
by definition. By using additionally the relation for the Christoffel symbols, see formula 
( |2.11[ ), we get 

bj = ^ ] tljidnj = y ^ bfi Qni 9 ~ ~ ^ ^ '^n l^m ^ ^ 9 ' dnniT — ~ ^ ^ ^nm l^m 

n n,i n,m i n,m 

for all J = 1, . . . , d. With v^{^^, oj) = kj this shows 

on TU . And due to uj ■ Vui = J2j '^j djT^s) ■ V^j = J2j '^j gfj we conclude 

on TU. Next we write 

gradM$(0 =(l- J[ {J^ J[y' J^) = E djT 

j 

for some uniquely determined Cj(^) G M. Hence it holds J2n'^n9in = • diT for each 
z = 1, . . . , d, since Sir = 0. Consequently, 

Cj = E '^nj = E 5ni 9'' = E ' ^iT = J2 9^' || 

for all J = 1, . . . , d. Altogether, we obtain 

^ ^ dq^ dvi 

on TU . Now we choose some r £ {!,..., A^}. We write 

limmer) = (l- J[ {J^ J[y^ Jc) er = E ^1 d^r 

i 

for some uniquely determined dl = d[{S,,uj) gM, i = 1, . . . ,d. With a similiar calculation 
as above we obtain d[ = J2j 9*"' ■ djT for all i = 1, . . . , d. Thus 

Mr = a(l- J[ ( J^y^ J^ e, • = a E d[ dir-V^ = aY, 9'' • djT ^. 
Consequently, we get 

92 



^ / ^ ^ \ ' ' 7 /euc\ ' ' 'euc Qyi-Qy'n- 

ij,n,m 



This finally implies 



QyiQyn z-^ dv'^dvi 
and the desired formula is proven. □ 
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Remark 5.2. In order to simulate the geometric Langevin process numerically, we may 
of course consider SDE (4.22) in M^'^ directly. Nevertheless, this requires a numerical 

So let us describe an alternative simulation method: Note 
as above, we may directly write down a 



algorithm which stays on TP 
that for a given coordinate system Tq of TP 
stochastic differential equation in the obvious way for each such local coordinates having 



the (local) Langevin operator L from Proposition 5.1 as associated Kolmogorov operator. 
In other words, this so called local SDE gives a (local) L-diffusion, thus yields a convenient 
way to simulate the (weakly unique) geometric Langevin diffusion process obtained from 
(4.22). We illustrate this method in the following applications. 



6. The Spherical Langevin process 

In practice there are many kinetic models availabe dealing with manifold-valued Brow- 
nian motions, especially spherical ones. Nevertheless, for many applications this gives to 
rough paths and it is reasonable to search for smoother versions of the classical models. 
Therefore, let us consider the following general point of view: The geometric Langevin 
process on TM, or more precisely its ^-coordinates, may be seen for general potential <I> 
as some kind of a smooth analogon to the Ornstein-Uhlenbeck process on M. The latter 
is given by the Stratonovch SDE on M of the form 

^2 



d6 



gradM^'(Ct) dt o-nM[Ct] ° dW*. 



(6.27) 



Here W denotes a standard A^-dimensional Brownian motion. In particular, for <I> = 0, 
Equation (6.27) gives a Brownian motion on M (with diffusion constant a), for details see 
e.g. |Hsu02| Prop. 3.2.6.] and |IIsu02| Theo. 3.1.4.]. In this situation, the corresponding 
Langevin process on TM with $ = can be used as the smooth analogon to the Brownian 
motion on M and yields the development of new, extended kinetic models, see for example 
the smooth fiber lay-down model developed in Section [7| In the latter section we are 
dealing with examples involving the spherical situation, i.e., M = S^. In view of these 
applications and in order to illustrate the strategy described in Remark 5.2 we shall 
seperately discuss now the geometric Langevin process living on TS'^. 

Therefore, consider the special coordinate charts {U, q) of given as the inverse of 
T : {a,a + 27r) x (0, vr) M, a e M, where 



(cos ^1 sin ^2 J sin 6*1 sin 02i cos ^2) , ^1 G (a, a -|- 27r), 6'2G(0,7r). 



By using formula (2.11) one easily calculates 



gn 



12 



sm 



21 



522 - 
cot 62 , 



11 



,22 



sm' 



11 



sm &2 cos 6/2, 



and else. So q^ = 9i, = 62 and = ki, = K2- Thus a local (Ito or Stratonovich) 
SDE simulating the local L-diffusion constructed from (4.22) reads as follows (for abuse 
of notation we omit the time index t) : 



d^i = Ki dt, d02 = dt 

dni = —Xki dt — 2 cot 92K1K2 dt — 

dK2 = —Xk2 dt -|- sin 62 cos 62^1 dt 



(6.28) 



1 

sin2 e2d9i' 



'dt + a- 



sm c'2 



dWl 



(1) 



802 



'1,C2 



idt + a dW, 



(2) 
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Figure 1. Spherical Brownian motion 



Here W denotes a standard two-dimensional Brownian motion and due to the 2tt- 
periodicity of r in the variable 6i, the state space of the previous SDE may be considered 
as M/27rZ x (0,7r) x M^. In this case the SDE shall more precisely be understood as an 
abstract Stratonovich equation similiar as in Section |4.1[ The latter can be simulated via 

directly. 

Via 



considering Equation (6.28) in M x (0,7r) x 

But first, we do a transformation of ( |6.28[ ) on R/2ttZ x (0,7r) x ™2 



i\, U2, Ki sm t/2, K2) 



With these new local coordinates, an easy calculation shows that (6.28) can equivalently 
be formulated on M/27rZ x (0,7r) x ^ 



d9i 
di/2 



sm 02 
— dt 



dt, (162 = 1^2 dt 

cot 92l'll'2 dt 



(6.29) 



1 (9$ 



6'2)dt + (jdW, 



(1) 



-Xi'2 dt + cot ^21^1 dt 



sin 62 dOi 

|J(01,02)dt + (7dW^i'\ 



These new local coordinates exactly are obtained by the (more natural) parametrization 
of TS^ given through 



where ni 



^/27rZ X (0,7r) x 9 (6*1,6*2,2^1,1^2) ^ {r,uini + 1^2^12) 

rdg^T and 722 = dg^r are the spherical unit vectors. Compare this with 



the general parametrization from (2.10). Morever, a stationary solution to the (formal) 



Fokker-Planck equation associated with the local SDE (6.29) now takes the desired well- 
known form e~ e ^2 ig ^2 2gjj^^2. 

A local SDE for the spherical Brownian motion with diffusion constant a is obtained 
via an SDE on M/27rZ x (0, vr) of the form 

,2 



d^i 



a 



(1) 



d02 



y cot 62 dt + adwf\ 



sm ^2 

This is due to the fact that the Kolmogorov operator associated with the latter equation 
coincides with the generator ^A§2 from ( |6.27 ) (for <I> = 0) restricted to ]R/27rZ x (0,7r), 
see again jHsu02| Theo. 3.1.4]. Figure [T] shows the paths of spherical Brownian motion 
with (T = 1. Compare the rough paths with the smooth ones from Figure [2] where the 
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Figure 2. Spherical Langevin process 

.^-coordinates of the spherical Langevin process are plotted for $ = 0, A = 1 and a = 1. All 
simulations, now and in the following, are always obtained with a simple Euler-Maruyama 
scheme. 

7. A NEW SMOOTH FIBER LAY-DOWN MODEL INVOLVING THE SPHERICAL LANGEVIN 
PROCESS AND FURTHER APPLICATIONS TO BIOLOGY 

In this section we present a new smooth three-dimensional fiber lay-down model. It 
is fascinating to note that the stochastic part of this model is given by the spherical 
Langevin process living on TS^ introduced previously. And moreover, as decribed in the 
introduction, we even present possible applications in biology to self-propelled interacting 
systems describing the collective behaviour of swarms of animals. 

7.1. A new smooth fiber lay-down modeL First remember the basic fiber lay-down 
Stratonovich SDE (|3.16l) with state space R"' x S"'"^ For convenience, we stay general 



and let d E N, (i > 2. The aim is to develop a fiber lay-down model having smoother 
trajectories than the original one, see the motivation in Section [3] Thus basically, we 
have to replace the spherical Brownian motion (with diffusion constant a) therein through 
the spherical Langevin process with potential equal to (and with diffusion constant a) 



from (4.21). Nevertheless, the latter lives on TS'^ ^. Consequently, the state space of the 
smooth fiber lay-down model must be equal to x TS'^^"'^. In notation 

R'^ X TS'^^^ = {(^, oj, £ R'^ X R'^ X R"^ \ = 1, u- fj, = 0}. 

Thus the vector field ( ^ nw^/^nY iC,^) & R"^ x S'^'^, with $ G C°^(M^) only 

V-(/-a;(g)w)V^>(0/ \ J J 

depending on ^, from the basic model has to be extended in the most natural way to some 

vector field ^ on M'^ x TS'^~^. We propose the ansatz 



A{C, uj,fx) := \-{I-u}(g> w)V$(0 , (C, OJ, fi) eR"" X TS 

G(e,a;,/i) ; 



.1 



UJ 

d . , qpfcifi— 1 



Now by condition (|2]7|, or see also ( |4.19[ ), the unknown G G C°°(M'' x TS'^"^) has to be 



chosen in order that 
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This means 

G{C, oj,fJ-)-uj = V^>(0 -/i for OJ, ji) gR'^x TS'^-i which is satisfied by setting 
G{C,uJ,fi) := /X • V^>(0 u for each £ R'^ x TS'^-^. So the smooth fiber lay-down 

model is defined with help of the following Stratonovich SDE as 

d^ = tjdt (7.30) 

doj = -{I - OJ (j) V^>(^) dt + /U dt 

d/i = /X • V^>(,^) wdt - A/xdt - |/x|^a;dt + cr (/ - cj (g) tj) o dWf 

with state space M"' x TS"'"^. Here is a standard d-dimensional Brownian motion and 
A, cr are nonnegative constants. For notational convenience, the time index t is omitted 
in the SDE. As already mentioned in the introduction, the model has been developed in 
cooperation with the authors from |KMW12b] . It has first been detected in its full form 
including all necessary terms by the second named author of the underlying article. 

Remark 7.1. Alternatively, one can proceed as follows. An equivalent form of the two- 



(7.31) 



dimensional basic fiber lay-down equation (3.16) is given by 

d^ = t(0) dt 

d6 = -T^{e) ■ V$(0 dt + dWt 

with W being a one-dimensional standard Brownian motion, t{6) := {cos{9) , sm{9))^ , 
T"*- := 1^ and $ S C°°{R?). The equation has more precisely to be understood as an 
abstract Stratonovich SDE with state space R? x M/27rZ and equivalent means here that 
any solution coincides weakly with any solution to ( 3.16| ) (for d = 2), see |GKMSl"2] for 
details. Now a natural smoother version of (7.31 ) reads 

dC = t{9) dt (7.32) 

de = -T^(0) • V$(Odt + z^dt 

du = -Xi^dt + adWt 

with state space x M/27rZ x R and A some nonnegative constant. To derive a natural 
higher dimensional version of the latter one proceeds exactly as in the Langevin case, see 
Section |4| Therefore, consider the diffeomorphism 

X R/2ttZ X R 9 (^, 9, v) ^ (C, OJ, n) 




with image x TS^. Now by following identically the argumentation from Section 4.2 
and by additionaly calculating the pushforward vector field of 

(/ - u}(g)uj)V^\ 



-T ■ simply as 



fl ■ V$ LJ J 



we obtain that the natural higher dimensional version of (7.32) is 



given by (7.30). In this sense both mathematical derivations are consistent and commu- 
tative. 



described in Remark 15.2 



Nevertheless, Equation (7.30) is not yet in suitable form for numerical simulations. As 



we need to derive a so called local SDE. For the original basic 



fiber lay-down equation (3.16) this is done in [ GKMS12] in each dimension. Therefore, we 



shall also calculate such equations for the smooth model from above for each d G N, d > 2. 
Moreover, this then also gives a geometrically precise justification of the calculations done 
in |KMW12b] and contains the practical relevant smooth three-dimensional fiber lay-down 



model in local form as special case. The proof partly follows the one from Proposition 5.1 



20 



MARTIN GROTHAUS AND PATRIK STILGENBAUER 



Therefore, we introduce first the foUowing special local parametrization of S", n G N, 
given by : M/27rZ x (0,7r)"-i ^ S" with t^^\9i) := (cos(ei), sin(0i))^, Oi G M/27rZ, 
and inductively 



r(")(ei,...0„) 



V"-i)(0i,...,0„_i)sin(0„) 
cos(6'„) 



, n G N, n > 2. 



In the following we write r instead of r^"). Note that = for i j. Furthermore, the 
spherical unit vectors on S", n G N, are defined by 

1 dr r~- dr 



re,- 



\deM de, 



do 



, j = l,...,re. 



(7.33) 



For the rest of this section, let Y := M*^ x M/27rZ x (0, vr) '^"^ x M^"^ with d as before and 
Y = X M/27rZ x M in case d = 2. Then a local parametrization of x TS'^-i is given 

by 



Y9(e,^,z^)^(e,u;,/x) := 




Finally, define 



A 



Ig^'^'de^rin ■ re,- 



for i,n,j = 1, . . . J d — 1. We arrive at the desired result. Recall our definition of a local 
SDE, see Remark [S^j and let L : C°°(M'^ x TS^^^^) C°°(M'^ x TS"^-^) be the Kolmogorov 
operator associated to the smooth fiber lay-down model. 



Proposition 7.2. A local SDE on Y associated to the smooth fiber lay-down SDE (7.30) 
generating also a (local) L- diffusion process is given by 



d^ = T{e) dt 

d0j = -^fg^ V^>(C) • nj{9) dt + Uj dt 

dvj = ^ /\inj V^>(0 • ni{e) dt - ^ /\inj Vi i^n dt - A i/j- dt + £7 ml^^ 



(7.34) 



where j = 1, . . . , d — 1. is a standard [d — 1)- dimensional Brownian motion and 
$ G C°°(M'^). /re particular, for d = 2 Equation (7.34) reduces to (7.32) and for d = ?i we 
obtain the three-dimensional local SDE for the smooth fiber lay-down model. 

Proof. At first, we get the pushforward vector fields 





d_ 



d_ 



where is the i-th unit vector in W^. Now there are uniqely determined a^, 6j, Cj G M, all 
depending on (^,a;,z^), such that 



d 



d 



_d_ 
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Figure 3. Smooth fiber- lay down 



Clearly Oj = for all i = 1, . . . ,d. And since (/ 



a;)V^>(0 = J2 j V^>(C) • rij nj we 



obtain 



I gjj V$(^) • rij for all j = 1, . . . , d — 1. By taking the scalar product with 



respect to some ^ on both sides of the last equation this yields 



i 



h ^ fn de^Un - rij = Y^ Aj„j V$(^) • rii Un- 

n 



Next, the previous arguments can be repeated with redefined ai,bj,Cj for the smooth 
vector field from 



X TS*^ ^ given by 




uj, n) gR'^ X TS' 



d-i 



Then Oj = for all i = 1, . . . ,d, and bj 



'g^i Uj for all j = 1, . . . , d — 1. Hence 



„ Aj„j Vi Vn- Note that the remaining vector fields in the smooth fiber lay down- 



model are already computed in the proof of Proposition 5.1 Following the notation from 
the latter, we have 

d ^ d 



This implies 



d_ 



Altogether, the generator associated to the local SDE from the proposition coincides with 
the generator of the smooth fiber lay-down model computed on Y. □ 

In Figure [3] we plot examplified the ^-trajectory of the two- and three-dimensional 
smooth fiber lay down model via simulating Equation (7.34). Therefore, for d = 3 one 
easily calculates A112 = — cot02) A121 = cot ^2 and else. The simulations are obtained 
by setting A = 1, ^> = and cr = 0, (T = 0.5, cr = 3. Finally, we remark that in analogy 
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to the basic fiber lay-down model, the stationary distribution of the associated Fokker- 
Planck equation corresponding to the smooth fiber lay-down model is again explicitly 
known, see |KMW12b] . Parameters as well as the potential can then be chosen such 
that the stationary distribution coincides with the real fiber-distribution as observed in a 
realistic industrial fiber lay-down process. 

Remark 7.3. We do not aim here to discuss the physical relevant scenarios for the choice 
of o" and A in case of a realistic fiber lay-down process as observed in the production 
process of nonwovens. Instead, for the smooth model we refer to |KMW12b] . See also 
|KMW12a] as well as |KMW09] and references therein. Moreover, we remark that an 
additional parameter B can be included into the three-dimensional model for measuring 
the anisotropic orientation of the fibers. This can be done as presented in |KMW12a] and 
|KMW12b] . An industrial relevant discussion of the existing (basic and smooth) fiber lay- 
down models is planned in a forthcoming research article of the authors from |KMW12b] 
together with the authors of the underlying paper. 

7.2. Self-propelled interacting particle systems with roosting. In |CKMTl"0] a 

two-dimensional spherical velocity model about self-propelled interacting particle systems 
with roosting force describing the collective behavior of swarms of animals is derived. The 
equation reads (at first informally) as 

The index i means the i-th particle with i = 1, . . . , K, K £ N and W = (Wi, . . . , Wk) is a 
i^-dimensional standard Brownian motion. Here r and T"*" are defined as above, r S (0, oo) 
is the constant velocity, a S (0,oo) the diffusion constant, $ and V are suitable potential 
functions. For details, derivation and motivation we refer to (CKMTlOj . Via defining 

^iiCi) ■=m) + ^T.i^{\^j-^i\), i = i,...,K, 

the two-dimensional model reads 

dCi = rT{0i)dt (7.35) 

d9i = -l T^i9i) ■ Ve,^(eO dt + ^ dWi^t 

where i = 1, . . . , K . So the model has an analogous form as the original two-dimensional 
fiber lay-down model and should precisely be understood as an abstract Stratonovich 
SDE with state space (M^ x R/2ttZ)^ , see |GKMS12] for the differential geometric details. 
Now by following again [GK MS12| . the higher dimensional (coordinate free) version of 
this model has state space (M'^ x S^)^, d G N, (i > 2, and is simply given by 

dCi = oji dt (7.36) 

duji = -UsM] V^^'fiCi) dt + Us^[uJi\ o dWi,t 

where i = 1, . . . , K. Here W = {Wi, . . . , Wk)'^ is a d • ii'-dimensional standard Brownian 



motion and S,- denotes the sphere of radius r > in M . In case d = 2, Equation (7.36) 
reduces to the two-dimensional model (7.35), see |GKMSl"2] for details (and compare this 



also with the first example in Section 9|) . For d = 3 one gets a physical relevant three- 
dimensional model, see |GKMSl"2] for an explicit numerical simulation formula. Moreover, 
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in exactly the same way as in our smooth fiber lay-down model from Section 7.1 



one 



obtains a smoother version of (7.36) now in the obvious way. Its applicability for physical 



relevant situations has of course to be checked and is left for future research. 



8. The geometric Langevin equation with spherical velocity 



Next we aim to derive a natural analogue of the geometric Langevin process (4.22) 
moving with velocity of constant absolute value, i.e., |a;| = r, r > 0. This is motivated by 
the mathematical structure detected in the fiber lay-down model in Section [3| In order to 
derive this Langevin model with velocity of constant absolute value, we follow the strategy 
as described in the fiber lay-down case in Section |3] and project the coordinates uj in the 
geometric Langevin equation (4.22) onto the sphere of radius r given by {w | |a;p = r^}. 
Then the resulting process should have state space S^M. Here recall that §,.M denotes 
the spherical tangent bundle bundle of M of radius r. In the following, we realize this 
equation in mathematical precise form. 

Therefore, note first that the orthogonal projection from onto the submanifold 



{OJ € 



pN 



r } is given by 



1 



jv-i a; 



I UJ. 



w e S, 



N-l 



We abbreviate ^ by and for r = 1 we shortly write S instead of Si. By the discussion 
above, it is natural to propose the following Stratonovich SDE 



ujf dt 



^.37) 



F{Ct,uJt) - gradM $(6)) dt + a (UsM] Hm^]) ° dWt. 



with some A^-dimensional standard Brownian motion W and state space S^.M, r > 0. 
Of course, this equation must be understood more precisely as an abstract Stratonovich 



equation similiar as the geometric Langevin equation, see (4.22). In the equation above 



ns^[a;]nM[C] denotes the usual matrix product between n§,.[ti;] and nM[C]- is 
defined in (4.24) and ^ E C°°(M). One easily verifies that the vector fields involved in 
(8.37) are really tangential to S^-M. Indeed, observe that given a vector field A satisfying 



condition (2.7), then the vector field B with 



fulfills conditions (2.7) and (2.8). Furthermore, we have ns,,[a;] (Aw) 

OJ 



\uj]J^X = jjx — {uJ., Jjx^ = Jjx — ( UJ, x) 



as well as 

= X 



for all {i,uj) G S^M and each x G M'^. Hence we have n§^[6t;]F(^,6t;) = F{^,uj) for each 
(■^,1^) G S^M. Moreover, recall the definition of gradM<I>(^) from (2.13). Altogether, the 
analogue manifold-valued Stratonovich SDE for the geometric Langevin process moving 
with velocity of constant absolute value simply reduces to 

dit = ujt dt (8.38) 
duJt = uJt) dt - (n§Ja;t]nM[6]) V$(6) dt + a (HsJw^] Uml^t]) o dWt 



38) the geometric Langevin equation with 



and has state space S^M. We call Equation 

spherical velocity. Note that in case M = M*^ it comes out the fiber lay-down Stratonovich 
SDE (3.16), see Section [sj Here recall the identification made in (2.14). Further examples 
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of this process are presented in the last section in order to convince the reader from this 
beauty geometric nature of the stochastic kinetic equation and to give further apphcations. 

Remark 8.1. It is fascinating to note that this equation arised from industrial real world 
applications. Indeed, the second named author of this article got this idea via trying to 
understand the geometry occuring behind the fiber lay-down world. 

9. Examples of the geometric Langevin process with spherical velocity 

In this section we provide some concret examples of our geometric Langevin equation 
moving with velocity of constant absolute value. At the one hand, this visualizes the dy- 
namics and shows its wonderful geometric nature. At the other hand, the second example 
below can be used in specific applications. Therein, we discuss the geometric Langevin pro- 
cess with spherical velocity on S^M where M = S^. Analogously to the spherical Langevin 
process discussed in Section [S] this stochastic process (or more precisely its ^-coordinates) 
serves as a further smooth process on S^. In particular, it can also be used to extend 
stochastic kinetic models having to rough paths due to an appearing Brownian motion 
on S^. A concrete example of this is again our basic (three-dimensional) fiber lay-down 
model, see Section [3j Consequently, a further smooth three-dimensional fiber lay-down 
model can be constructed. As mentioned in the introduction, we intend to publish an 
additional research article concerning this application. 

9.1. The cylinder M = x M. We start with the cyhnder M = x M. Thus 

pS TipS I /-2 I t2 _ -I ( il\ ( _ n I, ,|2 _ J2^ 



S,,M = {(e,^) G X Uf + e2 = 1, (^J^j • = 0, = r'], r > 0. 

In this section let Y := M/27rZ. There exists a natural diffeomorphism between Y x M x Y 
and SrM given by 



/ V \ 

rri (f- 
V rT2 ) 



Here _L denotes the first derivative w.r.t. Q or a respectively. For abuse of notation we 
omit sometimes the arguments of the occuring functions as well as the time index t. We 
have the following result. The proof works similiar as the one from the upcoming result 
in our second example, see the Appendix of this section. 

Lemma 9.1. The generator L : C°°(S,.M) C°^{8rM) with r > of the geometric 
Langevin equation with spherical velocity on S^M, M = x M, is given on Y x M x Y as 

d a 1 , _ d 

Here ^(a,^) := $(^(a), ^s), (a, 6) e Y x M /or $ G C 



'OO / 



Consequently, an equivalent SDE on Y x R x Y generating the (weakly unique) L- 
diffusion process reads 

d(a, ^3)^ = rT(0) dt (9.40) 

de = - T^iO) ■ V^(a, 6) dt + - dWi . 
r r 
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Figure 4. Langevin process with spherical velocity on S x M 



Remark 9.2. This result is not surprising, having the explicit two-dimensional fiber lay- 
down equation in mind, see e.g. Section[7]or |GKMSl"2] . Assuming a periodic boundary on 
the latter in the first variable, we obtain ( |9.40 ). This even shows that we can alternatively 
observe the general structure of the spherical velocity version of the geometric Langevin 
equation at our original fiber lay-down equation with the help of the previously defined 
state space transformation (9.39 1 . 

In Figure |4] we plot the ^-coordinates of the geometric Langevin equation with spherical 
velocity on SM (r = 1) via simulating (9.40). We choose <I> = and cr = on the left 
hand side and a = 4 on the right hand side. 

9.2. The sphere M = S^. Now we set M := S^. Then 

S,,M = {(^, cj) G X I = 1, ^ • w = 0, |a;|2 = r^}, r > 0. 
A local parametrization of S^M is given on Y x (0, vr) x Y, Y := M/27rZ, as 



Yx(O,7r)xY9(0i,e2,a)^(e,w) = ( 
In this example r and ip are defined as 



ripi ni + rip2 n2 



t{0i,(^2) '■= (cos ^1 sin 6*2, sin 6*1 sin 02, cos 6*2)'^ , ^{a) := (cos a, sin a)'^ 

and ni = ni{9i), 712 = 71.2(^1, 6*2) are the spherical unit vectors ni := \de-^T\'^^ Oq-^t , n2 := 
de^T, see Equation ( 7.33[ ). We have the following result, for the proof we refer to the 
appendix of this section. 

Lemma 9.3. The generator L : C°°(SrM) C°°(SrM) with r > of the spherical velocity 
version of the geometric Langevin equation with state space S^M, M = S^, is given locally 
on Y X (0, vr) X Y as 



cos a d d 
sm^2 001 062 

^ 1 sina „ ^ 

r cos a cot 62 -\ tt-t Oq, w 

r sin(6'2) 



(9.41) 



cos a do^ ^' 



d 



+ 



a 



52 



da 2r2 da"^ 



In this case ^(^1,^2) := ^(r(6'i, ^2)), (^1,^2) e Y x (0,7r), where $ E C°°(J 
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Figure 5. The geometric Langevin process with spherical velocity on 



In particular, an SDE modelling the (local) L-diffusion on Y x (0, vr) x Y reads 



COSQ , 

dOi = r — dt 

sm 62 

A62 = r sin a dt 



(9.42) 



da = r cos a cot 6*2 dt H — — Bq, dt cos a 80^^ di -\ — dWt ■ 

r sm O2 T f 

Remark 9.4. Note that a stationary solution to the (formal) Fokker-Planck equation 



associated to (9.42) is given on Y x (0,7r) x Y by e ^'^ sin02- This form is expected 
having the well-known form of the stationary solution to the fiber-lay down Fokker-Planck 
equation in mind, see Section [3} 

In Figure |5] we plot the ^-coordinates of the geometric Langevin equation with spherical 
velocity on SM (r = 1) via simulating Equation (9.42). a increases from cj = 0, a = 0.1, 
(J = 1 to cr = 4. Always <^ = 0. 



9.3. Appendix. We shall give the proof of Lemma |9.3| since the spherical Langevin pro- 
cess with spherical velocity can be of interest in specific applications, see the motivation 
from the introduction of this section. 
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Proof of Lemma 9.3. As before A denotes the pushforward of some vector field A. We 
have 

^ ^ / sin 6*2 ni ^ A = "2 A ^ ^ / 

Let X G be arbitrary. Now there exists uniquely determined oi, . . . , 03 G M with 

\ ^ ^ ^ 



nsja;]nMK]x; -"^a^i +"'9^2 ^°'aa- 

Immediately, we conclude that ai and a2 must be zero. By taking the scalar product with 
respect to on both sides of the last equation we get 

03 = ^ (n§,,MnM[^]a;, -(/?2?^l + '/'l"'2)euc 

= ^ (nM[e]a^, -V2ni + y'in2)euc = ^ ('^^' (x • ng) 
This follows since (a;, —ip2ni + 93ln2)g^^. = and nM[^]a; = x ■ ni ni + x ■ n2 112. Here 



09 = (—(f2,fi) ■ Let = r , r^, , ej the i-th unit vector in M'^, i = 1, . . . ,3. 
By the previous calculation Ai writes on Y x (0, vr) x Y as 



d 



Ai = - yip ,Zij — , Zi:= \ I , ? = 1, . . . ,3. 



r ^ ^ euc oa \ei ■ n2 

Consequently, for each such i we have 



1 / \ / ± \ d 1 / ± \2 5^ 

^ / euc ^ / euc Oa ^ ^ cue Oa'^ 



An easy calculation shows X)?, ('^"'"i-^j') = 1. Hence fv''''' -^^i) = 0- The latter 

^ ' euc ^ ' euc 

means [^iZi] ((/?-'", Zj) = 0. So altogether X^i-^? = ^W^- Next observe that 

^ ^ euc ^ ^ euc 

— ^' = sin 02^1 -V^^, —^ = n2-Ve.^. 
So the vector field Vi := — \ „ r ^-J^ r^n„ , ) writes on Y x (0,7r) x Y in the form 

^ \ de.,^ ) da 

For the remaining vector field we again make the ansatz 



^^^=1 FliM}^''W,^^'W2^^'L 
for some uniquely determined hi G M, i = 1,...,3. Here recall the definition of F in 
Section jsj Thus directly bi = and 62 = ''^<^2- Furthermore, one easily checks that 

F{^,uj) = r^^ = r. By taking again the scalar product with respect to on both 
sides of the last equation we obtain 

= 61 {ifi dg^ni + cos 92ip2 ni, -ip2ni + ipin2)g^^ + 63 = -h cos 62 [ifl + + 63- 
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And therefore 63 = 61 cos ^2 = r cos a cot 6*2- Altogether, the claim follows since the the 



generator L is equal to L = Vi + V2 + ^ Ya=i -^f • 



□ 
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